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Abstract. Studies of single spin asymmetries for pion electroproduction in semi- inclusive deep- 
inelastic scattering are presented using the polarized ~6 GeV electrons from at the Thomas 
Jefferson National Accelerator Facility (JLab) and the Continuous Electron Beam Accelerator 
Facility (CEBAF) Large Acceptance Spectrometer (CLAS) with the Inner Calorimeter. The cross 
section versus the azimuthal angle (ph of the produced neutral pion has a substantial sin ())h amplitude. 
The dependence of this amplitude on Bjorken xb and on the pion transverse momentum is extracted 
and compared with published data. 



1. Introduction 

In recent years it became clear that understanding the orbital motion of partons, and in particular 
the role of partonic initial and final state interactions of quarks, is crucial for the construction of a 
more complete picture of the nucleon in terms of elementary quarks and gluons. Parton distribution 
functions have been generalized to contain information not only on the longitudinal momentum 
but also on the transverse momentum distributions of partons in a fast moving hadron. Intense 
investigation of Transverse Momentum Dependent (TMD) distributions of partons both from the 
experimental and theoretical sides, indicate that QCD-dynamics inside hadrons is much richer than 
what can be learned from collinear parton distributions. 

Two fundamental mechanisms have been identified leading to single-spin asymmetries (SSAs) in 
hard processes: the Sivers mechanism pi [21 [3l IU |5] , which generates an azimuthal asymmetry in 
the distribution of quarks in the nucleon due to their orbital motion, and the Collins mechanism 
[HE], which generates an asymmetry during the hadronization of quarks. TMD distributions also 
contain unique information on the role of initial and final state interactions of active partons in 
hard scattering processes O HI [5] . 

Semi-inclusive deep-inelastic scattering (SIDIS) has emerged as a powerful tool to probe nucleon 
structure and provide direct access to TMDs through measurements of spin and azimuthal 
asymmetries. In particular, SSAs, arising from correlations of the transverse parton momentum and 
the transverse spin of the parton or from the initial or final state hadron, provide unprecedented 
information about spin-orbit correlations. QCD factorization for semi-inclusive deep-inelastic 
scattering at low transverse momentum in the current-fragmentation region has been established 
in references [3 |8l [9]. This new framework provides a rigorous basis to study the TMD 
parton distributions from SIDIS data using different spin-dependent and independent observables. 
The analyses of the TMDs strongly depends also on knowledge of the fragmentation functions 
[10^ [TTl [T2[ [13] . Many experiments are currently trying to pin down various TMD effects through 
1) semi-inclusive deep-inelastic scattering ( HERMES at DESY [H [El [H [IT], COMPASS at 



CERN \T8], CLAS and Hall- A at Jefferson Lab [13 EOl [2T]), 2) polarized proton-proton collisions 
(PHENIX, STAR and BRAHMS at RHIC [22l [231 [M]) and 3) electron-positron annihilation (Belle 
at KEK [25J). 

The cross section for single pion production by longitudinally polarized electrons scatteried from 
unpolarized protons may be written in terms of a set of structure functions [6l [26l I27j . The helicity 
dependent part {(Tlu) arises from the anti-symmetric part of the hadronic tensor: 



The subscripts on ctlu specify the beam and target polarizations, respectively (L stands for 
longitudinally polarized and U for unpolarized). The azimuthal angle (j)h is the angle between 
leptonic and hadronic plane according to the Trento convention [28]. The kinematic variables 
xb, y, and z are defined as: xb = Q'^ /2{Piq), y = {Piq) / {Piki) , z = {PiP)/{Piq), where 

= —q^ = —{ki — k2f' is the four-momentum of the virtual photon, ki {k2) is the four-momentum 
of the incoming (scattered) lepton, Pi and P are the four momenta of the target nucleon and the 
observed final-state hadron (respectively), 7^ = AM'^ x'^y'^ / Q'^ , M is the nucleon mass, and Ae is 
the electron beam helicity. 

The structure function F^jj'^''' arises due to the interference of the longitudinal and transverse 
photon contributions. In the partonic description of SIDIS, assuming factorization holds, 
contributions to structure functions can be written as convolutions of parton distribution and 
fragmentation functions dependent on the scaling variables xb and z, respectively [6l[7]. 

The beam-spin asymmetries in single-pion production off an unpolarized target are higher-twist 
by nature |29l 130] . Higher-twist observables are important for understanding the long-range quark- 
gluon dynamics and may also be accessible as leading contributions through the measurements of 
certain asymmetries [261 [311 [32l |33] , in particular the beam SSAs. 

Recently, higher-twist effects in SIDIS were interpreted in terms of average transverse forces 
acting on the active quark at the instant after being struck by the virtual photon [34]. Different 
contributions to the beam SSAs discussed so far provide information on leading and sub-leading 
parton distribution and fragmentation functions, related both to Collins and Sivers production 
mechanisms. Sizable beam SSAs were predicted [30j based on spin-orbit correlations as the 
dynamical origin. Within this framework, the asymmetry generated at the distribution level, is 
given by either the convolution of the T-odd parton distribution with the twist-3 fragmentation 
function E [31], or the convolution of the twist-3 T-odd distribution function with the unpolarized 
fragmentation function Di[35]. 

2. The TT*^ beam spin asymmetry 

In this section we present measurements from the EOl-113 CLAS dataset of beam-spin asymmetries 
in the electroproduction of neutral pions in deep-inelastic scattering using the 5.776 GeV electron 
beam and the CEBAF Large Acceptance Spectrometer (CLAS) [36] at Jefferson Laboratory. 
Longitudinally polarized electrons are scattered off a liquid- hydrogen target. The beam polarization, 
frequently measured with a M0ller polarimeter, was on average 0.80 with a fractional systematic 
uncertainty of 3%. The beam helicity was flipped every 30 ms to minimize systematic instrumental 
effects. Scattered electrons were detected in CLAS. Electron candidates were selected by a 
hardware trigger using a coincidence between the gas Cherenkov counters and the lead-scintillator 
electromagnetic calorimeters (EC). 

Neutral pion events were identified by calculating the invariant mass of two photons detected 
with the CLAS electromagnetic calorimeter and the Inner Calorimeter (IC) [37]. In each kinematic 
bin, TT*^ events are selected by a gaussian plus linear fit (see figure [Tj). The combinatorial background 
is subtracted in each bin from the number of events inside 2>a using the linear component of the fit. 
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Figure 1. Invariant mass spectrum of the 
77 system M(77) in an arbitrarly chosen xb, 
Prp, z and (/)/i-bin, fitted by gausian plus hnear 
polinomial. Vertical black lines indicates ±3a" 
from the mean. 



Figure 2. Example of a po " sm{(j)h) ht for 
0.4 < 2 < 0.7, 0.1 < XB < 0.2 and 0.2 
GeV/c< Pt < 0.4 GeV/c. Only statistical 
error bars are shown. 



Deep-inelastic scattering events are selected by requiring > 1 GeV^ and W'^ > 4 GeV^ , where 
W is the invariant mass of the hadronic final state. Events with low missing mass of the en^ system 
{Mx < 1-5 GeV) were discarded to exclude contributions from exclusive processes. The minimum 
value of the transverse vr'' momentum, Pt > 0.05 GeV, ensures that the azimuthal angle (j)h is well 
defined. The total number of selected electron- vr*^ coincidences is ~ 3.0 x 10^. 

The beam-spin asymmetry has been calculated for each kinematic bin as: 



P N^,{<^h) + 

where P = 0.794 ± 0.024 is the absolute beam polarization this data set and N\-) and A'^ o sxa 
the number of vr'^'s with positive and negative beam helicity, respectively. Asymmetry moments 
are extracted by fitting the i;^ft-distribution of Alu in each xb and Pt bin with the theoretically 
motivated function pQ ■ sm((j)h). An example of this fit is shown in figure [2] for an arbitrarly chosen 
kinematic bin. 

In figure [3l the extracted moment is presented as a function of Pt for different xb ranges. 

Systematic uncertainties, presented by the bands, include the uncertainties due to the background 
subtraction, the event selection and possible contributions of higher harmonics in the extraction of 
the moments. These contributions are added in quadrature. An additional 3% uncertainty due to 
the beam polarization measurement and another 3% uncertainty from radiative effects should be 
added to the systematic uncertainties. 

The ^^f/*^'" moment increases at low Pt and reaches a plateau at values of about 0.3 GeV. There 
is an indication that the decrease of ^j^^ at large Pt-, expected from perturbative QCD, starts 
already at Pt ~ 0.6 GeV. 

The surprising characteristics of favored and unfavored Collins functions, being of roughly equal 
magnitude but having opposite signs, as indicated by latest measurements at HERMES [TU [IB] . 
COMPASS [IB] and Belle [25], puts the vr*^ in a unique position in SSA studies. Since the vr^ 
fragmentation function (FF) is the sum of vr"'" and vr" EEs, its favored and unfavored contributions 
will be roughly equal and, in the case of the Collins FF, will cancel each other to a large extent. 
Contributions to the beam-SSA related to spin-orbit correlations can thus be studied without a 
significant background from the Collins mechanism. 
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Figure 3. ^'/■j/ as function of Pt for different xb ranges and integrated over 0.4 < z < 0.7. 
The error bars correspond to statistical and the bands to systematic uncertainties. An additional 
3% uncertainty arises from the beam polarization measurement and another 3% uncertainty from 
radiative effects. 



The measured beam-spin asymmetry amplitude for vr^ appears to be comparable with the vr^ 
asymmetry from a former CLAS data set [38j both in magnitude and sign, as shown in figure [H 
indicating that contributions from the Collins mechanism cannot be the dominant ones. 

A similar measurement has been performed by the HERMES collaboration at a higher beam 
energy of 27.6 GeV [17]. After taking into account the kinematic factors in the expression of the 
beam-helicity dependent and unpolarized terms ([27]) 



j^y> = 1 — , 2/0 ' ^) 

CLAS and HERMES measurements are found to be consistent as shown in figures [5] and O indicating 
that at energies as low as 4-6 GeV the behavior of beam-spin asymmetries is similar to higher 
energy measurements. The CLAS data provide significant improvements in precision of beam SSA 
measurements in the kinematic region where the two data sets overlap, and extend the measurements 
to the large xb region not accessible by HERMES. 



2.1. Contamination from two-hadron production 

Li this section a study of the beam spin asymmetry of single-pion production from exclusive and 
non-exclusive vector meson production is presented. In order to investigate a possible contamination 
of the presented results by pions originating from the decay of exclusively produced vector mesons. 




Figure 4. The vr*^ beam-spin asymmetry amplitude ^^Jy"^'" as function of xb compared to that for 
7r+ from an earlier CLAS measurement [38] . Uncertainties of the measurement are as in figure [3l 
For both data sets < Pt >w 0.38 and 0.4 < z < 0.7. 
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Figure 5. ^^Ic/ multiplied by the kinematic 
factor < Q > / f{y) as a function of xb from 
CLAS and HERMES [IT]. The 0.4 GeV/c< 
Pt < 0.6 GeV/c range for the CLAS data was 
used to compare with HERMES, since this is 
the closest kinematic range. 



Figure 6. multiplied by the kinematic 

factor < Q > / f{y) as function of Pt from 
CLAS and HERMES [17J (the same as in 
figure [5]). The 0.1 < xb < 0.2 range for the 
CLAS data was used to compare with, since 
this is the closest kinematic range. 



as advocated in jl7], events with at least one vr*^ and ir'^ were selected from the same dataset. 
Exclusivity was ensured by requiring the missing mass Mx of the ep — t- e'iT^7r~^X system to be 
within 0.8 GeV < Mx{ep — ?■ e'7r^7:~^X) < 1.1 GeV. The invariant mass distribution of this exclusive 
two pion sample is presented in the upper panel of figure [7] where the exclusive peak is clearly 
visible along with non-resonant exclusive two-pion production. The vertical lines indicate invariant 
mass ranges for which the beam spin asymmetry amplitude A^^^'^'' was extracted for each of the two 
exclusively produced pions following the extraction method as described before. These amplitudes 
are presented for vr'^ and vr^ in the second top panel of figure [7] as labeled therein, along with the 
average values of < Pt > and < z > for the invariant mass ranges. 

Figure [8] is the same as figure [7| except for two-hadron events with Mx{ep — >• e'lT^TT^x) > 1.5 
GeV (non exclusive). 

The ^ijj''^ values for 7r'^(7r'^) increase (decrease) with increas invariant mass, and reach their 
maximum absolute value at the highest invariant mass. 




Figure 7. From top to bottom: invariant 

mass M^07r+ of two hadrons, ^^Ic/'^'' ' ^^"^ 
average < Pt > and < z > as function of 
M^o^+ for TT^ and 7r+ in exclusive two-pion 
production. 



Figure 8. From top to bottom: Invariant 
mass M^o^+ of two hadrons, A^^^'^'' , and 
average < Pt > and z as function of M^Ot^+ 
for vr'^ and vr^ in non-exclusive production 
{M^{ep e'TT^TT+x) > 1.5 GeV). 



In contrast to a similar study for presented in |17j , the asymmetres for vr'' and vr"*" have minimal 
values in the invariant mass range where events dominant, and furthermore the highest bin of 
invariant mass, the asymmetries for vr" and tt^ have different signs. It should be noted also, that 
in these studies the momentum distributions and the acceptance for the two pions are not similar, 
while in [17] the two pions have similar momentum distributions and acceptance. ^^Ic/"^'' °^ 
SIDIS ir^ and tt"*" in the same kinematic bin have the same sign and their amplitudes are roughly 
equal (figured]). For the case of two-hadron production (figure [8]), in all invariant mass bins (except 
the first) A^^'^''^ for vr*^ and 7r+ have different signs, which is another reason for careful studies of 
two hadron production. 

Comparison of the measured beam SSAs for pions originating from electroproduction of two-pions 
in exclusive and non-exclusive processes, suggests that the SSAs depend mainly on the kinematics 
of the pion, in particular on z and Pt, and are comparable for pions from different samples within 
the same kinematical bin. 

3. Conclusions 

We have presented measurements of kinematic dependences of the beam-spin asymmetry in semi- 
inclusive tt" electroproduction from the EOl-113 CLAS dataset. The sin (ph amplitude is extracted as 
a function of xb and the transverse pion momentum Pt, integrating over the z-range 0.4 < z < 0.7. 

The asymmetry shows no significant x^-dependence over the measured range. The Pt 
dependence of ^^Ic/"^'' consistent with an increase at low values of Pt which reaches a plateau for 



values Pt > 0.3 GeV. There is an indication that the decrease of A^^ij at large Pt, expected from 
perturbative QCD, already starts at Pt ~ 0.6 GeV. The observed asymmetry amplitudes for vr^ 
indicate that the major contribution pion beam SSAs may be due to spin-orbit correlations. 

The results obtained are compared with published HERMES data [T7], providing significant 
improvement in precision and an important input for studies of higher twist effects. Despite the fact 
that the partonic formalism is much better suited for higher energy reactions, there is a reasonable 
agreement in size and behavior of beam SSAs measured over a wide energy range. 

Preliminary results on two-hadron production show interesting trends: opposite signs of 
asymmetries for vr'^ and tt"*", the highest asymmetry for the highest invariant mass range. The 
beam spin asymmetry from vr*^ and vr"*" two-hadron production has been studied for the first time 
as a function of the invariant mass of two pions. The A^li/''^ amplitude of vr*^ from exclusive and 
non-exclusive two-pion production exhibits a similar magnitude and similar kinematic dependences 
as for semi-inclusive tt^ production. Therefore, a possible contribution of this exclusive channel to 
our semi-inclusive result would not alter the presented asymmetry amplitudes. 
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